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Abstract —Wireless device-to-device (D2D) communication un¬ 
derlaying cellular network is a promising concept to improve user 
experience and resource utilization. Unlike traditional D2D com¬ 
munication where two mobile devices in the proximity establish a 
direct local link bypassing the base station, in this work we focus 
on relay-aided D2D communication. Relay-aided transmission 
could enhance the performance of D2D communication when 
D2D user equipments (UEs) are far apart from each other 
and/or the quality of D2D link is not good enough for direct 
communication. Considering the uncertainties in wireless links, 
we model and analyze the performance of a relay-aided D2D 
communication network, where the relay nodes serve both the 
cellular and D2D users. In particular, we formulate the radio 
resource allocation problem in a two-hop network to guarantee 
the data rate of the UEs while protecting other receiving nodes 
from interference. Utilizing time sharing strategy, we provide a 
centralized solution under bounded channel uncertainty. With a 
view to reducing the computational burden at relay nodes, we 
propose a distributed solution approach using stable matching 
to allocate radio resources in an efficient and computationally 
inexpensive way. Numerical results show that the performance 
of the proposed method is close to the centralized optimal solution 
and there is a distance margin beyond which relaying of D2D 
traffic improves network performance. 

Index Terms —Device-to-device (D2D) communication, LTE- 
A L3 relay, uncertain channel state information, distributed 
resource allocation, stable matching. 

I. Introduction 

We consider relay-assisted device-to-device (D2D) commu¬ 
nication underlaying LTE-A cellular networks where D2D user 
equipments (UEs) are served by the relay nodes HI, 12- When 
the link condition between two D2D UEs is too poor for 
direct communication, the D2D traffic can be transmitted via a 
relay node, which performs scheduling and resource allocation 
for the D2D UEs. We refer to this as relay-aided D2D 
communication which can be an efficient approach to provide 
better quality of service (QoS) for communication between 
distant D2D UEs. For this, we utilize the LTE-A Layer-3 
(L3) relay 0. We consider scenarios in which the potential 
D2D UEs are located in the same macrocell; however, the 
proximity and link condition may not be favorable for direct 
communication. Therefore, they communicate via relays. The 
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radio resources (e.g., resource blocks [RBs] and transmission 
power) at the relays are shared among the D2D communication 
links and the two-hop cellular links. 

We formulate an optimization problem to allocate radio 
resources at a relay node in a muti-relay multi-user orthogonal 
frequency-division multiple access (OFDMA) cellular network 
(e.g., LTE-A network). Due to the NP-hardness of the re¬ 
source allocation problem, we utilize time sharing strategy 
and provide an asymptotically optimal centralized solution. 
Considering the random nature of wireless channels, we re¬ 
formulate the resource allocation problem using the worst-case 
robust optimization theory. The uncertainties in link gains are 
modeled using ellipsoidal uncertainty sets. Each relay node 
can centrally solve the problem taking channel uncertainty into 
consideration. However, considering the high (e.g., cubic to 
the number of UEs and RBs) computational overhead at the 
relay nodes, we provide a distributed solution based on stable 
matching theory which is computationally inexpensive (e.g., 
linear with the number of UEs and RBs). We also analyze the 
stability, uniqueness, and optimality of the proposed solution. 

Considering the computational and signalling overheads 
and lack of scalability of the centralized solutions, game 
theoretical models have been widely used for wireless resource 
allocation problems. However, the analytical tractability of 
equilibrium in such game-theoretical models requires special 
properties for the objective functions, such as convexity, which 
may not be satisfied for many practical cases 10. In this 
context, resource allocation using matching theory has several 
beneficial properties 0, 0- For example, the stable matching 
algorithm terminates for every given preference profile. The 
outcome of matching provides suitable solutions in terms of 
stability and optimality, which can accurately reflect different 
system objectives. Besides, with suitable data structures, a 
Pareto optimal stable matching (e.g., allocation of resources to 
the UEs) can be obtained quickly for online implementation. 

The goal of this work is to design a practical radio resource 
scheme for relay-aided D2D communication in a multi-user 
multi-relay OFDMA network. As opposed to most of the work 
in the literature where channel gain information is assumed 
to be perfect, we capture the dynamics of random and time- 
varying nature of wireless channels. To this end, we develop 
a low-complexity distributed solution based on the theory 
of stable matching and demonstrate how this scheme can 
be implemented in a practical LTE-A system. The major 
contributions of this work can be summarized as follows: 

• We model and analyze the radio resource allocation 
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problem for relay-aided D2D communication underlaying 
an OFDMA cellular network considering uncertainties in 
channel gains. We formulate an optimization problem to 
maximize system capacity in a two-hop network while 
satisfying the minimum data rate requirement for each 
UE and limiting the interference to other receiving nodes. 
We show that the convexity of the optimization problem 
is conserved under bounded channel uncertainty in both 
the useful and interference links. 

• Using the theory of stable matching, we develop a 
distributed iterative solution, which considers bounded 
channel uncertainty. The stability, uniqueness, optimality, 
and complexity of the proposed solution are analyzed. 
We also present a possible implementation approach of 
our proposed scheme in an LTE-A system. 

• Numerical results show that the proposed distributed solu¬ 
tion performs close to the upper bound of the optimal so¬ 
lution obtained in a centralized manner; however, it incurs 
a lower (e.g., linear compared to cubic) computational 
complexity. Through simulations, we also compare the 
performance of the proposed approach with a traditional 
underlay D2D communication scheme and observe that 
after a distant margin, relaying of D2D traffic improves 
network performance. 

We organize the rest of the paper as follows. We briefly 
review the related work in Section [TT] Section [III] presents 
the system model and related assumptions. Followed by the 
formulation of the nominal resource allocation problem in 
Section [TV] we reformulate the resource allocation problem 
considering wireless channel gain uncertainties in Section [V] 
We develop the stable matching-based distributed resource 
allocation algorithm in Section [VI] Theoretical analysis of 
the proposed solution is presented in Section IV III In Section 
IVIIII we present the performance evaluation results before we 
conclude the paper in Section [IX] 

II. Related Work 

Despite the fact that the resource allocation problems in 
cellular D2D communication have been intensively studied 
in the recent literature, only a very few work consider re¬ 
lays for D2D communication and incorporate wireless link 
uncertainties in the formulation of the resource allocation 
problem. In (6), a greedy heuristic-based resource allocation 
scheme is proposed for both uplink and downlink scenarios 
where a D2D pair shares the same resources with cellular 
UE (CUE) only if the achieved signal-to-interference-plus- 
noise ratio (SINR) is greater than a given SINR requirement. 
A resource allocation scheme based on column generation 
method is proposed in 0 to maximize the spectrum uti¬ 
lization by finding the minimum transmission length (i.e., 
time slots) for D2D links while protecting the cellular users 
from interference and guaranteeing QoS. A new spectrum 
sharing protocol for D2D communication overlaying a cellular 
network is proposed in ED. which allows the D2D users 
to communicate bi-directionally while assisting the two-way 
communications between the eNB and the CUE. A graph- 
based resource allocation method for cellular networks with 


underlay D2D communication is proposed in 0. A two- 
phase resource allocation scheme for cellular network with 
underlaying D2D communication is proposed in iflOl . In I'm, 
the mode selection and resource allocation problem for D2D 
communication underlaying cellular networks is investigated 
and the solution is obtained by particle swarm optimization. 
The works above do not consider relays for D2D communi¬ 
cation. 

A distributed relay selection method for relay-assisted D2D 
communication system is proposed in El- In HD, d, 
authors investigate the maximum ergodic capacity and outage 
probability of cooperative relaying in relay-assisted D2D com¬ 
munication considering power constraints at the eNB. Taking 
the advantage of L3 relays supported by the 3GPP standard, 
in 0 , a centralized resource allocation approach is proposed 
for relay-assisted D2D communication assuming that perfect 
channel information is available. A gradient-based distributed 
resource allocation scheme is proposed in 0 for relay-aided 
D2D communication in a multi-relay network under uncertain 
channel information. In this multi-relay network, the interfer¬ 
ence link gain between a UE and other relays (to which the UE 
is not associated with) is modeled with ellipsoidal uncertainty 
sets. However, the uncertainty in direct channel gain between 
relay and the UE is not considered. In this paper, we remodel 
the previous formulation and extend the work in 0 by 
incorporating uncertainties in both the useful and interference 
links. In particular, we present a distributed resource allocation 
algorithm using stable matching considering the uncertainties 
in wireless channel gains (e.g., channel quality indicator [CQI] 
parameters according to the LTE-A terminology). 

Although not in the context of D2D communication, match¬ 
ing theory has been used in the literature to address the radio 
resource allocation problems in wireless networks. A spectrum 
allocation algorithm using matching theory is proposed in ini 
for a cognitive radio network (CRN) under perfect channel 
assumption. In El, a two-sided stable matching algorithm 
is applied for adaptive multi-user scheduling in an LTE-A 
network. A distributed matching algorithm is proposed in El 
for cooperative spectrum sharing among multiple primary and 
secondary users with incomplete information in a CRN. In 
El, a distributed algorithm is proposed to solve the user 
association problem in the downlink of small cell networks 
(SCNs). A matching-based subcarrier allocation approach is 
proposed in |[T9l for services with coupled uplink and down¬ 
link QoS requirements. The radio resource (e.g., subcarrier 
and power) allocation problem for a full-duplex OFDMA net¬ 
work is modeled as a transmitter-receiver-subcarrier matching 
problem in l20l . 

The matching-based solutions proposed in m-m do 
not consider D2D-enabled networks. In the context of D2D 
communication, most of the works (e.g., 0, 0, 0, 0- 
031, El, El) provide centralized solutions. Also, note that 
in 0-ED, the effect of relaying on D2D communication is 
not investigated. Moreover, the wireless link uncertainty is 
not considered in 0, 0-El, El-Go). Different from the 
above works, we propose a stable matching-based distributed 
radio resource allocation approach considering the channel 
gain uncertainties in a multi-relay and multi-user relay-aided 




3 


TABLE I 

Summary of Related Work and Comparison with Proposed Scheme 


Reference 

Problem focus 

Relay-aided 

Channel 

information 

Solution approach 

Solution type 

Optimality 

m 

Resource allocation 

No 

Perfect 

Proposed greedy heuristic 

Centralized 

Suboptimal 

m 

Resource allocation 

No 

Perfect 

Column generation based 
greedy heuristic 

Centralized 

Suboptimal 

0 

Resource allocation 

No* 

Perfect 

Numerical optimization 

Semi-distributed 

Pareto optimal 

m 

Resource allocation 

No 

Perfect 

Interference graph coloring 

Centralized 

Suboptimal 

0 a (M 

Resource allocation 

No 

Perfect 

Two-phase heuristic 

Centralized 

Suboptimal 

° -B rn 
a § 

Resource allocation, 
mode selection 

No 

Perfect 

Particle swarm optimization 

Centralized 

Suboptimal 

§ § 03 

■g 1 

Resource allocation, 
mode selection 

Yes 

Perfect 

Proposed heuristic 

Distributed 

Suboptimal 

11 m 

Theoretical analysis, 
performance evaluation 

Yes 

Perfect 

Statistical analysis 

Centralized 

Optimal 

0 

Performance evaluation 

Yes 

Perfect 

Heuristic, simulation 

Centralized 

N/A* 

m 

Resource allocation 

Yes 

Perfect 

Numerical optimization 

Centralized 

Asymptotically 

optimal 

iD 

Resource allocation 

Yes 

Uncertain* 

Gradient-based iterative update 

Distributed 

Suboptimal 

03 

Resource allocation 

N/aT 

Perfect 

One-to-one matching 

Centralized 

Optimal 


Cross layer scheduling 

N/A+ 

Perfect 

Many-to-one matching 

Centralized 

N/A* 

Jj O 

3 OB 

Spectrum sharing 

N/At 

Complete, 

incomplete 

One-to-one matching 

Distributed 

Pareto optimal 

m-B (7b) 

Cell association 

N/At 

Perfect 

Many-to-one matching 

Distributed 

N/At 

O 

Resource allocation 

N/At 

Perfect 

Many-to-one matching 

Distributed 

N/A* 

£ HQ 

Resource allocation 

N/At 

Perfect 

One-to-one matching 

Centralized 

N/A* 

Proposed 

scheme 

Resource allocation 

Yes 

Uncertain 

Matching theory (many-to-one 
matching) 

Distributed 

Weak Pareto 
optimal 


tNo information is available. 

^^Uncertainty in channel gain in the direct link between UEs (relays) and relays (eNB) is not considered, 
t Not applicable for the considered system model. 


D2D communication scenario. A summary of the related work 
and comparison with our proposed approach is presented in 
Table U 

III. System Model 

A. Network Model 

Let C = {1, 2, ..., L} denote the set of fixed-location L3 
relays in the network (Fig. 1 in [0). The system bandwidth is 
divided into N orthogonal RBs denoted by A f = {1,2,..., N} 
which are used by all the relays in a spectrum underlay 
fashion. The set of CUEs and D2D pairs are denoted by 
C = {1,2,..., C} and D = {1,2,..., D}, respectively. We 
assume that association of the UEs (both cellular and D2D) 
to the corresponding relays are performed before resource 
allocation. Prior to resource allocation, D2D pairs are also 
discovered and the D2D session is setup by transmitting known 
synchronization or reference signals EH . 

We assume that the CUEs are outside the coverage region of 
the eNB and/or having bad channel condition, and therefore, 
the CUE-eNB communications need to be supported by the 
relays. Communication between two D2D UEs requires the 
assistance of a relay node due to poor propagation condition. 
The UEs assisted by relay l are denoted by Ui. The set of 
UEs assisted by relay l is Ui = {1,2 ,,U{\ such that Ui C 
{C U £>}, \/l £ £, U; Ui = {C U Dj, and f), Ui = 0. 

In the second hop, there could be multiple relays transmit¬ 
ting to their associated D2D UEs. We assume that multiple 
relays transmit to the eNB (in order to forward CUEs’ traffic) 
using orthogonal channels and this scheduling of relays is done 


by the eNB0. In our system model, taking advantage of the 
capabilities of L3 relays, scheduling and resource allocation 
for the UEs is performed in the relay nodes to reduce the 
computational load at the eNB. 


B. Achievable Data Rate 

(n) 

Let 7 Ui i j denote the unit power SINR for the link between 
UE ui £ Ui and relay l using RB n in the first hop and 7 ^ 2 
be the unit power SINR for the second hop. Note that, in 
the second hop, when the relays transmit CUEs’ traffic (i.e., 
ui £ {C l~l Ui}), 2 denotes the unit power SINR for the 

link between relay l and the eNB. On the other hand, when 
a relay transmits to a D2D UE (i.e., ui £ {DC Ui}), <y|; a 
refers to the unit power SINR for the link between relay l and 
the receiving D2D UE for the D2D-pair. 

(n) 

Let > 0 denote the transmit power in the link 

between i and j over RB n and Brjs is the bandwidth of 
an RB. The achievable data rat^B for ui in the first hop 
can be expressed as = B rb log 2 (l + 

Similarly, the achievable data rate in the second hop is 
r i" 2 = b RB log 2 (l + ■ Since we consider a two- 

1 Scheduling of relay nodes by the eNB is not within the scope of this work. 

2 We will present the rate expressions in Section flV-AI 
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hop communication, the end-to-end data ratif] for ui on RB 
n is half of the minimum achievable data rate over two hops 
, i.e., 

Ri.'SS}- <*> 


R { u? = 2 mm • 


r ( n ) - 


1, if RB n is assigned to UE ui 
0 , otherwise. 


( 2 ) 


N 0 denotes thermal noise. The interference link gain between 
relay (UE) i and UE (relay) j over RB n in hop k is denoted 
by g) j k , where UE (relay) j is not associated with relay (UE) 
i. Similarly, the unit power SINR for the second hofQ. 


IV. Resource Allocation: Formulation of the 
Nominal Problem 

In the following, we present the formulation of the resource 
allocation problem assuming that perfect channel gain infor¬ 
mation is available. This formulation is referred to as the 
nominal problem since the uncertainties in channel gains are 
not considered. 

For each relay, the objective of radio resource (i.e., RB 
and transmit power) allocation is to obtain the assignment of 
RB and power level to the UEs that maximizes the system 
capacity, which is defined as the minimum achievable data rate 
over two hops. Let the maximum allowable transmit power for 
UE (relay) is P™ ax (I y [ nax ) and let the QoS (i.e., data rate) 
requirement for UE ui be denoted by Q UI . The RB allocation 
indicator is a binary decision variable x'u? G { 0 , 1 }, where 


(n) 

%UI, 2 = 


E 

Vti,e{ DDUj}, 


n l,Ui,2 

<")p(») 


9j^NB,2 + cr ' 


■, uj € {c ni/J 


I.W 

‘."I . 2 


E 

iAUec 


1 3,U^»3,«J,2 


Ui G {TnWi} 


(5) 

where hi >Ult 2 denotes the channel gain between relay-eNB link 
for CUEs (e.g., ui G {C f) Ui}) or the channel gain between 
relay and receiving D2D UEs (e.g., ui G {V C\Ui}). From 
0 , the maximum data rate for UE ui over RB n is achieved 
when S ,1 = P& Therefore, in the second hop, 
the power Pi Ul allocated for UE «,/, can be expressed as 
a function of power allocated for transmission in the first 

hop, P Uul as follows: P/Jj = 

>l,u l ,2 

Hence, the data rate for ui over RB n can be expressed as 
= \B rb log 2 (l + ^ 1 ) ■ Considering the above, 

the objective function in 0 can be rewritten as 


'LjTi p O) 

WTZ “<>*• 


A. Objective Function 

N 

Let R Ul = x denote the achievable sum-rate 

n —1 

over allocated RB(s). We consider that the same RB(s) will 
be used by the relay in both the hops (i.e., for communication 
between relay and eNB and between relay and D2D UEs). 
The objective of resource allocation problem is to maximize 
the end-to-end rate for each relay l G £ as follows: 


max 


N 

w 

p(») 

i' i,u t uiEUi n=l 


E 


) 

Ul - LL Ul 


(3) 


N 

E + (6) 


“1 ’ ui £Ui n—1 


For each relay l £ £ in the network, the 

objective of resource allocation problem is to 


obtain 

the 

RB and 

power 

allocation 

xj 

= 


T (l) 

*^1 5 * ‘ 

r m 

• r^l j • 

r (i) 

■ ’ X Ui > • ■ 

Pi = 

A" 


P w 

l r l,l ’ 

p(l) 

• * ^U t ,l 

P W 
• ■ ■ ' r u u i 

which maximize 

the data rate. 



and 


r Wl T 
> x u, 

T 

, respectively. 


where the rate of UE m over RB n 

SpBlog^l + P^E)' 
Brb log 2 (l + Pifyiwfl) t 
In 0, the unit power SINR for the first hop. 


R { u? = 2 min ' 


(n) 

Tm, 1,1 




v ^ n ) n^ 

r Un,i^u 

Vuj £Uj, 

jj£i,je£ 


ULj ,j^Uj ,1,1 


where h^ k denotes the direct link gain between node i and 
j over RB n for hop k G {1,2}, er 2 = NqB rb in which 


(4) 


B. Constraint Sets 

In order to ensure the required data rate for the UEs while 
protecting all receiver nodes from harmful interference, we 
define the following set of constraints. 

• The constraint in 0 ensures that each RB is assigned to 
only one UE, i.e., 

^xW<l, VnGAC (7) 

uiEUi 


• The following constraints limit the transmit power in each 


3 In a conventional D2D communication approach where two D2D UEs 
communicate directly without a relay, the achievable data rate for D2D UE 
u G T> over RB n can be expressed as R^ = Brb log 2 (l + V 


where 


h£l 




is the channel gain in the link 


vjeUu 

between the D2D UEs and U u denotes the set of UEs transmitting using the 
same RB(s) as u. 


4 According to LTE-A standard, the L3 relays are able to peform similar 
operation as an eNB. Besides, the relays in the network are interconnected 
through X2 interface for better interference management ED Since the relays 
can estimate the CQI values (and hence the interference level) using X2 
interface, it is straightforward to account for interference in 0 and 0 . 
Consequently, interference from other transmitter nodes (e.g., UEs associated 
to other relays in the first hop or other relays in the second hop) will appear 
as a constant term in 0 and 0. 
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of the hops to the maximum power budget: 

N 

T, x % )p £h p r x ’ yui&Ui (8) 

n— 1 
N 

E E*& )p S< p r x ■ m 

ui Gl4i n =1 

« Similar to ff24ll . we assume that there is a maximum 
tolerable interference threshold limit for each allocated 
RB. The constraints in (fTTTb and ( 1 1 I [ i limit the amount of 
interference introduced to the other relays and the receiv¬ 
ing D2D UEs in the first and second hop, respectively, to 
be less than some threshold, i.e.. 


V x (n) P (n) o (n) < / (n) Vn G TV 

ui&Ai 

V x (n) P, (rt) 0 (n) < / (n) \/n G Af 

/ , r l,ui*l,u \,2 — 2’ V "' ^ * /V ‘ 

ui EUi 


( 10 ) 

(ii) 


• The minimum data rate requirements for the CUE and 
D2D UEs is ensured by the following constraint: 


R Ul > Q Ul , g Ui . (12) 


• The binary decision variable on RB allocation and non¬ 
negativity condition of transmission power is defined by 

G {0,1}, P%j > 0, Vu, € Uu Vn € AT. (13) 

Note that in constraint (TITTb and (fill , we adopt the concept 
of reference user. For example, to allocate the power level 
considering the interference threshold in the first hop, each UE 
ui associated with relay node l obtains the reference user u* 
associated with the other relays and the corresponding channel 
gain g^i j for Vn according to the following equation: 

u* = argmax 6 Ui,j £ hj 6 C. (14) 

3 

Similarly, in the second hop, for each relay Z, the transmit 
power will be adjusted accordingly considering interference 
introduced to the receiving D2D UEs (associated with other 
relays) considering the corresponding channel gain gj u . 2 for 
Vn, where the reference user is obtained by 

u* = argmax g[ n ^ 2 , j f l, j € £, Uj G {D n Uj}. (15) 

Uj 


C. Centralized Solution 

Corollary 1. The objective function in © and the set of 
constraints in mm turn the optimization problem to a 
mixed-integer non-linear program (MINLP) with non-convex 
feasible set. Therefore, the formulation described in Section 
HE is NP-hard. 

A well-known approach to solve the above problem is to 
relax the constraint that an RB is used by only one UE by 
using the time-sharing factor ||25l . In particular, we relax the 
optimization problem by replacing the non-convex constraint 
x G {0,1} with the convex constraint 0 < xff) < 1. Thus 

(n) (n) — 

Xui represents the sharing factor where each x ui denotes the 
portion of time that RB n is assigned to UE ui and satisfies 


the constraint xffl < 1, Vn. Besides, we introduce a 

uiGUi 

new variable S < ' n> , = xfflP^ n ] > 0, which denotes the actual 
transmit power of UE ui on RB n |26| . Then the relaxed 
problem can be stated as follows: 


(P2) max V R u . 

Xu i uiGUi 

subject to 0. (fl2l > and 


N 

E 

n—1 


s£h < pt x ’ ^ 


N 


E \ 0 -( n ) o( n ) Dl 

uieUi n— 1 


E 

ui£Ui 


q{n) (n) < An) w 

S ui,l9 u *,LI — 1 th,V’ Vn 

n(") c( n ) „( n ) <- r(™) w„ 

H u l ,l t, u l ,l9l,u’l,2 < 1 th, 2. Vn 


ui&Ai 


(16) 

(17) 

( 18 ) 

(19) 

( 20 ) 
( 21 ) 


where 7^,1 = 


E s. 

VlLj^Uj, 


k \w 

u l U 

Cn) „(») 


, 1 +^ 


rr( n ) 

’’ H ui,l 


u i’ L . 

zu&~ 


and 


N 


r ui = y \x^b rb io g2 (1 


Lu i ~ / > 2 

n= 1 


— zptr~ 


Corollary 2. The objective function in ( 1761) is concave, the 
constraint in (E3 is convex, and the remaining constraints in 
0 , mm are affine. Therefore, the optimization problem 
P2 is convex. 


Since P2 is a non-linear convex problem, each relay can 
solve the optimization problem using standard algorithms such 
as interior point method l27l Chapter 11], The centralized 
optimization-based solution is summarized in Algorithm [l] 
Each relay locally solves the optimization problem P2 and 
informs the other relays the allocation vectors using X2 inter¬ 
face. The process is repeated until the data rate is maximized, 
e.g., Ri(t)-Ri(t- 1) < e for VZ, where R t (-) = ^ #«,(•) 

ui&Ai 

is the sum data rate for relay l obtained by solving the 
optimization problem at iteration (•) and e is a small value. 

The duality gap of any optimization problem satisfying the 
time-sharing condition becomes negligible as the number of 
RBs becomes significantly large. The optimization problem 
P2 satisfies the time-sharing condition, and therefore, the 
solution of the relaxed problem is asymptotically optimal ||28S . 

Given the parameters of other relays (e.g., Xj, Pj Vj f 
l,j G £), at each iteration of Algorithm [l] the allocation 
vectors (e.g., xj, P 1 ) obtained at each relay l provide locally 
optimal solution for l. In addition. Algorithm [l] allows the 
relays to perform allocation repeatedly with a view to finding 
the best possible allocation. If the data rate at the (t + l)-th 
iteration is not improved compared to that in the previous 
iteration t, the algorithm terminates, and the allocation at 
iteration t will be the resultant solution. Each of the iteration 
of Algorithm |T] outputs the solution of relaxed version of 
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Algorithm 1 Optimization-based resource allocation 
1 : Each relay l £ C estimates the CQI values from previous time 
slot and determines reference gains g^l\ 1 and g\ n (. 2 ,Vui,n. 

2: Initialize t := 0. 

3: repeat 

4: Update t := t + 1. 

5: Each relay l £ C: 

• solves the optimization problem P2 and multicasts the 
allocation variables X;, P; to all relay j ^ l, j £ C over 
X2 interface. 

• calculates the achievable data rate based on current 
allocation as Ri{t) := E R ^y 

u l CUi 

6: until data rate not maximized and t < T max . 

7: Allocate resources (i.e., RB and transmit power) to associated 
UEs for each relay. 


the original NP-hard optimization problem. Since the solution 
of relaxed problem gives us the upper bound, and at the 
termination of Algorithm [T] we obtain an upper bound of 
the achievable sum rate. 

V. Resource Allocation Under Channel 
Uncertainty 

For worst-case robust resource optimization problems, the 
channel gain is assumed to have a bounded uncertainty of 
unknown distribution. An ellipsoid is often used (e.g., |29| - 
ED) to approximate such an uncertainty region. 


A. Uncertainty Sets 

(n) 

Let the variable F ^ u l denote the normalized channel gain 
which is defined as follows: 


z 7 '( n ) A ^A 

Ul,Uj,l , (n) 7 

n ui,l,l 


\/uj 


( 22 ) 


(n) 

In addition, let F Vi \ denote the uncertainty set that describes 
the perturbation of link gains for ui over RB n. The normalized 
gain is then denoted by 


f(") = p(") i AF^ 

Ui,Uj,l Ui,Uj,l ' Ui,Uj,l 


(23) 


where F*) n ] , is the nominal value and A F^ n \ , is the 

'LL l . Uj • L 'LL l . 'LL j • L 

perturbation part. The uncertainty in the CQI values is modeled 
under an ellipsoidal approximation as follows: 




(n) 


Ul ,Uj , 


AF 


(n) 


Vuj^Uj, 

jAUec 

(24) 


A n ).\ 


where (> 0 is the uncertainty bound in each RB. Using 
we rewrite the rate expression for ui over RB n as 


R U? = \B RB 10g 2 ( 1 + — +a tn) 


(25) 


Wuj EUj , 


where ai?'* = 


^T 


and F' 


(n) 


Ul ,Uj ,1 


is given by 


B. Reformulation of the Optimization Problem Considering 
Channel Uncertainty 

Utilizing uncertainty sets similar to (l24t in the constraints 
(IT8li- (l20l i. the optimization problem P2 can be equivalently 
represented under channel uncertainty as follows: 


uiGtii 


P3) 

max 


min 


E 


(26) 


x£\s£ 
1 u i 

AF ( ” : 

U l ’ 


» 

’ ui£Ui 





A H {n) , 

u t ,l 


l,U. 

,2 



subject to 

0 , 03, 03 

1 , (121 li and 



N 








EEC 

rr( n ) 1 

<l 

) c(«) 

) °u, ,1 

< 

r>max 

r l 


(27) 

uiGUi n= 1 








E (4 

A* s ~ 

<1 

+ 

‘e' 

\ c( n ) 
J S Ul,l 

< 

r( n ) 

1 th,li 

Vn 

(28) 

ui&Ai 








' ff( n ) n ) 
^uijyituj ,2 ~ 

^ A < 


\ q( n ) 
) b u u l 

< 

A n ) 

1 th, 2 ’ 

Vn 

(29) 


E 


") |2 
l,Uj,l 1 

<0 


' (30) 


Vuj&Uj 

EE \ AH Stf <&i) 2 (3D 

ui£Ui n= 1 

E I A 4ti ! 2 ^ (^) 2 > Vn(32) 

ui&Ai 

E l A <W n El 2 < (d:;) 2 , Vn(33) 

ui&Ai 


where for any parameter y, y denotes the nominal value 
and Ay represents the corresponding deviation part; £ 2 1 , dE 

(n) U ‘ 

and £4 are the maximum deviations (e.g., uncertainty 
bounds) of corresponding entries in CQI values. In P3, R UI 
is given by 


N 

Rut = E \ X u, >B KI1 X 

n— 1 


log 2 1 + 


u i if 




( „) + .(«) 


I ■ (34) 


\/uj £Uj, 


The above optimization problem is subject to an infinite 
number of constraints with respect to the uncertainty sets and 
hence becomes a semi-infinite programming (SIP) problem 
' ll32l . In order to solve the SIP problem it is required to 
transform P3 into an equivalent problem with finite number 
> of constraints. Similar to HD, 1301, we apply the Cauchy- 
Schwarz inequality ll33l and transform the SIP problem. More 
specifically, utilizing Cauchy-Schwarz inequality, we obtain 
the following: 


V A F {n) ,S {n) -< 

/ J u l i u j J u j i3 


Vuj GUj , 


E i Ai tEi 2 E i<)i 2 


\ Vuj&Aj, 


Vuj GlAj , 
oAUec 


<&] 


\ 


E «!,-)• < 35 > 


Vuj GlAj , 
jjtijec 
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Similarly, 


N 




N , s 2 


£ £ 

Q (n) 

< 

^2/ A 

££«f 

(36) 

ui&Ai n= 1 



\ 

Ul 6 \lAl Tl —1 



q (n) 

S ui,l 

< 

e( n ) 

i/£(0 2 

(37) 

uiGUi 




V ui&Ai 


V A /^( n ) „( n ) 

/ j ui,l9l,u* ,2 

q(n) 

< 

d:; 

i/£ («) 2 - 

(38) 


■ i i i / ■ — 

ui€Ui V u i€Ui 


Note that, as presented in Section IV-A1 to tackle the un¬ 
certainty in channel gains, we have considered the worst-case 
approach, e.g., the estimation error is assumed to be bounded 
by a closed set (uncertainty set). Hence, from (l35l >- (l38l) . under 
the worst-case channel uncertainties, the optimization problem 
P3 can be rewritten as P4, where R Ul is given by (143b . The 
transformed problem is a second-order cone program (SOCP) 
ll27l Chapter 4] and the convexity of P4 is conserved as shown 
in the following proposition. 

Proposition 1. P4 is a convex optimization problem. 


are determined by the given preference profiles. That is, for 
both the RBs and the UEs, the lists of preferred matches 
over the opposite set are maintained. For each RB, the relay 
holds its preference list for the UEs. The matching outcome 
yields mutually beneficial assignments between RBs and UEs. 
Stability in matching implies that, with regard to their initial 
preferences, neither RBs nor UEs have an incentive to alter 
the allocation. 

A. Concept of Matching 

A matching (i.e., allocation) is given as an assignment of 
RBs to UEs forming the set of pairs (ui ,n) €Ui x A f. Note 
that a UE can be allocated more than one RB to satisfy its 
data rate requirement; however, according to the constraint in 
©, one RB can be assigned to only one UE. This scheme 
corresponds to a many-to-one matching in the theory of stable 
matching. More formally, we define the matching as follows 

urn 

Definition 1. A matching /// for V/ G C is defined as a 
function, i.e., pi : Hi U A/" —> Hi U J\f such that 


Proof: Using an argument similar to that in footnote[4] the 
objective function of P4 in ( 139b is concave. The constraints in 
©. (fTTl) . (fTi! are affine and the constraint in ( 1 1 2b is convex. In 
addition, the additional square root term in the left hand side of 
the constraints in (l40l >. ( |3 1 b . and (l42l ) is the linear norm of the 
vector of power variables S^\ with order 2 , which is convex 
ll27l Section 3.2.4]. Therefore, the optimization problem P4 
is convex. ■ 

P4 is solvable using standard centralized algorithms such 
as interior point method. The joint RB and power allocation 
can be performed similar to Algorithm [T] and an upper bound 
for the solution can be obtained under channel uncertainty. It 
is worth noting that solving the above SOCP using interior 

/ /_ _\ 3 ^ 

point method incurs a complexity of O 


(xf+Pi)' 


i) pfn) £UiC {0} and pi{n) G {0,1} 

ii) pi(tii ) G A/ - and pi{ui) G {1 , 2 ,..., k ui } 

where the integer k Ui < N, pi(ui) = n 3=> p{n) = uiforMn G 
A/”, Vu/ G Hi and Pj(-) denotes the cardinality of matching 
outcome Pj(-). 

The above definition implies that pi is a one-to-one match¬ 
ing if the input to the function is an RB. On the other hand, 
pi is a one-to-many function, i.e., pi(ui) is not unique if the 
input to the function is a UE. In order to satisfy the data 
rate requirement for each UE, we introduce the parameter n Ul 
denoting the number RB(s) which are sufficient to satisfy the 
minimum rate requirement Q UI . Consequently, the constraint 


at 


N 


each relay node where y denotes the length of vector y. 
Besides, the size of the optimization problem increases with 
the number of network nodes. Despite the fact that the solution 
from Algorithm [l] outputs the optimal data rate, considering 
very short scheduling period (e.g., 1 millisecond in LTE- 
A network), it may not be feasible to solve the resource 
allocation problem centrally in practical networks. Therefore, 
in the following, we provide a low-complexity distributed 
solution based on matching theory. That is, without solving the 
resource allocation problem in a centralized manner using any 
relaxation technique (e.g., time-sharing strategy as described 


in the preceding section), we apply the method of two-sided 



■ flW- 

stable many-to-one matching 041. 

H/ = 

d(1) 

L R u, 

0 m 

H u, J 


VI. Distributed Solution Approach Under 
Channel Uncertainty 

The resource allocation approach using stable matching 
involves multiple decision-making agents, i.e., the available 
RBs and the UEs; and the solutions (i.e., matching between UE 
and RB) are produced by individual actions of the agents. The 
actions, i.e., matching requests and confirmations or rejections 


in (0 is rewritten as = k U i , \/ui. Generally this 

n— 1 , , 

parameter is referred to as quota in the theory of matching [5J- 
Each user ui will be subject to an acceptance quota n ui over 
RB(s) within the range 1 < n Ul < N and allowed for matching 
to at most n Ul RB(s). The outcome of the matching determines 
the RB allocation vector at each relay l, e.g., pi = x/. 

B. Utility Matrix and Preference Profile 

Let us consider the utility matrix 11; under the worst-case 
uncertainty, which denotes the achievable data rate for the UEs 
in different RBs, defined as follows: 


(44) 


where 11 i[i,j] denotes the entry of i-th row and j-th column 
in 11/, and R is given by (l45l) . Each of the UEs and RBs 
holds a list of preferred matches where a preference relation 
can be defined as follows f36l Chapter 2], 


Definition 2. Let A be a binary relation on any arbitrary set 
3. The binary relation A is complete if for \/i,j G 3, either 
















X. 


u i ’ u t ,i uiEUi 


subject to 0, (fl2l >. (fI7] >, m and 


N 


E E<XV&. 

ui&Ai n —1 


N 


iV „ 

EE(«) sfl” 


uiEUi n =1 


E -(n) o(n) >(n) / \ ' 7 c(") V 

9uf,l,l 0 ui,l + S3 U , i 
ui&Ai V ujGWj 

E fr(n) (n) q (n) An) I f o(n)V 

n ui,l9l, u *,2 a ui,l + S4„ */ fWui.U 

ui&Ai V “i£Wi 


< 7^ Vn 


< 7^ Vn 

— 1 th.2' i v "- 


(40) 

(41) 

(42) 


AT 


/ 


= E o Xu? Brb log 2 


n=l 


AT 


sW 


1 + 


y' f( n ) o( n ) , A n ) 

Z-/ r Ul ,Uj ,l°Uj ,j ' >l u 

E (s%) 2 + 4?> 

Vu.-ew,, \ 


&i,i ec \ 

jAiJec 


= E o x ^ Brb log2 




1 + 




pzM q(n) , An) An) 
Z> yu i .l,l°u i ,j T n ui .l.l^lu. 

E «1,) 2 +^ 2 


V«3 G Uj , 


jjtt, jec 


(43) 


E l) = log 2 


/ 


V 


p( n )u( n ) 
r Ul ,1 l Ul ,1,1 


x^ A^i) An) p( n ) , u( n ) A n ) 
l-J dju 3 y Ui ,l,l r u U j ~T~ rb ui,l, 1S1 U . 

E (4 } ^) 2 + ^ 


Mu-j G Uj , 

dec \ 

jAUec 


(45) 


i cl j or j >i i or both. A binary relation is transitive if 
i f j and j >i k implies that i >i k for Vfc € E. 77ze binary 
relation f is a (weak) preference relation if it is complete and 
transitive. 

The preference profile of a UE ui £ Ui over the set 
of available RBs A f is defined as a vector of linear order 
(P Ul (N) = iXi[ui,i]i£tf. The UE ui prefers RB m to n 2 
if m Z ri 2 , and consequently, il;[u/,ni] > il/[u/,n 2 ]. 
Likewise, the preference profile of an RB n £ Af is given 
by Tniph) =ih\j,n] jeUl - 


C. Algorithm for Resource Allocation 

Based on the discussions in the previous section, we utilize 
an improved version of matching algorithm (adapted from |[37l 
Chapter 1.2]) to allocate the RBs. The allocation subroutine, 
as illustrated in Algorithm [2] executes as follows. While an 
RB n is unmatched (i.e., unallocated) and has a non-empty 
preference list, the RB is temporarily assigned to its first 
preference over UEs, i.e., u/. If the allocation does not exceed 
k ui , the allocation will persist. Otherwise, the worst preferred 
RB from uf s matching will be removed even though it was 


previously allocated. The iterations are repeated until there 
are unallocated pairs of RB and UE. The iterative process 
dynamically updates the preference lists and hence leads to a 
stable matching. 

Once the optimal RB allocation is obtained, the transmit 
power of the UEs on assigned RB(s) is obtained as follows. We 
couple the classical generalized distributed constrained power 
control scheme (GDCPC) f38l with an autonomous power con¬ 
trol method ll39l which considers the data rate requirements of 
UEs while protecting other receiving nodes from interference. 
More specifically, at each iteration t, the transmission power 
for each allocated RB is updated as follows: 



A(t — 1), if A(f - 1 ) < P^ r ' 

P ^\, otherwise 


(46) 


where 


p(n) m ' 

ui 


A(i 1) — 2 R Ul (t-i)_ 1 - r U |,2 




£ <> (s",+o,,) £ £ i 

\ Tl= i v 1 


.(») 

U l 


(47) 

(48) 
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Algorithm 2 RB allocation using stable matching 

Input: The preference profiles TV, (A f), tP n (Ui)\ \/ui € Ui,n € Af. 
Output: The RB allocation vector x;. 

1: Initialize X; := 0. 

tv , 

2: while (3u; with ^ a:«, < R UI ) or (3 n with a4, = 0, Vu; € 

n=l 

W; and tP n (W;) 7 ^ 0 ) do 

3: Um P '■ = most preferred UE from the profile ‘P n (iUi). 

4: Set X^u^p := 1. /* Temporarily allocate the RB */ 

N 

5: if E x ul P > s% P then 

3=1 

6: nip := least preferred resource allocated to u mp - 

7: Set := 0. /* Revoke allocation due to quota violation */ 

8 : end if 

N 

9: if E x ( p‘l p = n Ump then 

3=1 

10: nip := least preferred resource allocated to u mp . 

11: /* Update preference profiles */ 

12: for each successor flip of m p on profile (M) do 

13: remove hi p from T’u mp (A/’). 

14: remove u mp from JV ip (W;). 

15: end for 

16: end if 

17: end while 


and P:' 1 ] is obtained as 

ui 

Hi ;] = min (Hll > min (Hi )max , )) ■ ( 49 > 

In ( l49l >. the parameter Hi). ' s chosen arbitrarily within the 

, n , 5(n) , 5 (n)^* , (n) • . , 

range of 0 < and j, is given by 


( n ) 

™u,,l = mm 


r(") 




r(™) 
1 th,2 

T^T ’ ^(n) —(n) 






(«) 


(50) 


Based on the RB allocation, the relay informs the parameter 
and each UE updates its transmit power in a dis¬ 
tributed manner using (l46t . Each relay independently performs 
resource allocation and allocates resources to corresponding 
associated UEs. The joint RB and power allocation algorithm 
is given in Algorithm |3] 


Algorithm 3 Joint RB and power allocation algorithm 


Phase I: Initialization 

1 : Each relay l € C estimates the nominal CQI values from 


previous time slot and determines reference gains g 

ffS,2’ Vw i’ n - 


(n) 

u 7 * Al¬ 


and 


2 : Initialize t := 0, '■= —jfi — Vui,n and if based on CQI 

estimates. 

Phase II: Update 

3: for each relay l € £ do 

4: repeat 

5: Update t := t + 1. 

6 : Build the preference profile tP n (Wi) for each RB n € fif 

based on utility matrix and inform corresponding entries of 
it; to UEs. 

7: Each UE ui G Ui builds the preference profile V ui {j\f). 

8 : Obtain RB allocation vector using Algorithm^ 

9: Update the transmission power using d46l i for Vu;, n and 

update the utility matrix it;. 

10 : Inform the allocation variables x;, P; to each relay j 7 ^ 

l,j £ £ and calculate the achievable data rate based on 
current allocation as Ri(t) := E *«.(*)■ 


uieiti 

11 : until data rate not maximized and t < T max . 

12 : end for 
Phase III: Allocation 

13: For each relay, allocate resources (i.e., RB and transmit power) 
to the associated UEs. 


From Definition [3] [ 4 ] the matching //; is blocked by RB n 
and UE ui if n prefers ui to p.;(n) and either i) ui prefers n 
to some n € or if) /i;(u;) < n ui and n is acceptable to 

ui. Using the above definitions, the stability of matching can 
be defined as follows EH Chapter 5], 

Definition 5. A matching jit is stable if it is individually 
rational and there is no pair ( ui,n ) in the set of acceptable 
pairs such that ui prefers n to pi(ui) and n prefers Ui to 
pi(n), i.e., not blocked by any pair of agents. 

Proposition 2. The assignment performed in Algorithm [2] 
abides by the preferences of the UEs and RBs and it leads to 
a stable allocation. 

Proof: See Appendix [H ■ 


VII. Analysis of the Proposed Solution 

In the following, we analyze the performance of our pro¬ 
posed distributed resource allocation approach under bounded 
channel uncertainty. More specifically, we analyze the stability, 
optimality, and uniqueness of the solution, and its computa¬ 
tional complexity. 

A. Stability 

Definition 3. (a) The pair of UE and RB ( ui,n ) in Ui x J\f is 
acceptable if ui and n prefer each other (to be matched) 
to being remain unmatched. 

(b) A matching pi is called individually rational if no agent 
(i.e., UE or RB) j prefers to remain unmatched to p,(j). 

Definition 4. A matching /// is blocked by a pair of agents 
(i,j) if they each prefer each other to the matching they obtain 
by pi, i.e., i h Pi{j) and j P pi(i). 


Note that the allocation of RBs is stable at each iteration 
of Algorithm [3] Since after evaluation of the utility, the 
preference profile of UEs and RBs are updated and the routine 
for RB allocation is repeated, a stable allocation is obtained. 

B. Uniqueness 

Proposition 3. If there are sufficient number of RBs (i.e., N > 
Ui), and the preference lists of all UEs and RBs are determined 
by the Ui X N utility matrix il; whose entries are all different 
and obtained from given uncertainty bound, then there is a 
unique stable matching. 

Proof: See Appendix |B] ■ 

C. Optimality and Performance Bound 

Definition 6. A matching ji/ is weak Pareto optimal if there 
is no other matching pi that can achieve a better sum-rate, 
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i.e., fii(-) > /q(-), where the inequality is component-wise and 
strict for one user. 

Proposition 4. The proposed resource allocation algorithm is 
weak Pareto optimal under bounded channel uncertainty. 

Proof: See Appendix [C] ■ 

Corollary 3. Since xj* satisfies the binary constraint in 
(0, and the optimal allocation (xj*,Pj*) satisfies all the 
constraints in the optimization problem P4, for a sufficient 
number of available RBs, the data rate obtained by Algo¬ 
rithm |3] gives a lower bound of the solution under channel 
uncertainty. 

D. Convergence and Computational Complexity 

Proposition 5. The subroutine for RB allocation terminates 
after some finite number of steps T'. 

Proof: Let the finite set X represent all possible com¬ 
binations of UE-RB matching where each element xf 1 e X 
denotes that RB j is allocated to UE i. Since no UE rejects 
the same RB more than once (see line 7 in Algorithm [2]>, the 
finiteness of the set X ensures the termination of RB allocation 
subroutine in finite number of steps. ■ 

Note that the distributed approach replaces the optimiza¬ 
tion routine in the centralized approach with the matching 
algorithm (Algorithm 2). Since by Proposition [5] we show 
matching algorithm terminates after finite number of iterations, 
ultimately Algorithm 3 will end up with a (local) Pareto 
optimal solution after some finite number of iterations. 

In line 6-7 of Algorithm |3] the complexity to output the 
ordered set of preference profiles for the RBs using any stan¬ 
dard sorting algorithm is O (NUi log Ui) and for each UE, the 

complexity to build the preference profile is O (TV log TV). Let 

Ui _ 1 v _ 

n (Ui) = 2NUi be the total length 

ui — 1 n—1 

of input preferences in Algorithmic] where TP ; (-) denotes the 
length of the profile vector (■). From Proposition [5] and 
G7] Chapter 1] it can be shown that, if implemented with 
suitable data structures, the time complexity of RB allocation 
subroutine is linear in the size of input preference profiles, 
i.e., 0(/3) « O(NUi). Since Phase II of Algorithm |3] 
runs at most fixed T max iterations, at each relay_node l, the 
complexity of the proposed solution is linear in TV and Ui. 

E. Signalling Over Control Channels 

Assuming that the relays obtain the CQI prior to resource 
allocation, the centralized approach does not require any 
exchange of information between a relay node and the asso¬ 
ciated UEs to perform resource allocation. However, in the 
distributed approach, the relay node and the UEs need to 
exchange information to update the preference profiles and 
transmit power. In both the approaches, the relay nodes need 
to exchange the allocation variables among themselves (e.g., 
over X2 interface) in order to calculate the interference levels 
at the receiving nodes. 


In the distributed approach, the exchange of information 
between a UE and the relay node during execution of the 
resource allocation algorithm can be mapped onto the standard 
LTE-A scheduling control messages. For scheduling in LTE- 
A networks, the exchanges of messages over control channels 
are as follows ESI. The UEs will periodically sense the 
physical uplink control channel (PUCCH) by transmitting 
known sequences as sounding reference signals (SRS). When 
data is available for uplink transmission, the UE sends the 
scheduling request (SR) over PUCCH. The relay, in turn, 
uses the scheduling grant (SG) over physical downlink control 
channel (PDCCH) to allocate the appropriate RB(s) to the UE. 
Once the allocation of RB(s) is received, the UE regularly 
sends buffer status report (BSR) using PUCCH in order to up¬ 
date the resource requirement, and in response, the relay sends 
the acknowledgment (ACK) over the physical hybrid-ARQ 
indicator channel (PHICH). Given the above scenario, the UEs 
may provide the preference profile ‘J’ UI (A/") with the SR and 
BSR messages. The relays may provide the corresponding 

values in the utility matrix, e.g., u Ul i = iljiq, j]j=i ... n and 
- ( n ) max ’ 

inform the parameter P-fi using SG and ACK messages. 
Once the RB and power allocation is performed, the relays 
multicast the allocation information over X2 interface. 

In what follows, we analyze signalling overhead for our 
proposed solution. For a sufficient number of available RBs, 
we consider two cases: a) number of available RBs at each 
relay l is equal to the number of UEs (e.g., TV = Ui and 
k Ui = 1, V'«/); and b ) the number of RBs is greater than 
the number of UEs (e.g., TV > Ui). For the first case, once 
Algorithm 2 terminates, all RBs are allocated to the UEs. 
This is because, by the definition of individual rationality (see 
Definition 3), none of the agents (i.e., UE or RB) wants to 
remain unallocated. Hence, at the end of any iteration t of 
Algorithm 2, there are TV — t unallocated RBs at each relay l. 
Therefore, the maximum number of iterations, say T max are 
required when all the RBs are allocated, e.g., TV — T max = 0, 
and therefore, T max = TV. Since at each iteration t, TV — 
t + 1 messages are exchanged, the total number of messages 
exchanged in Algorithm 2 can be quantified as 

n, = T g;‘(jv-i + i) = JV(JV + 1) . (51) 

i—1 

In Algorithm 3, each relay l exchanges the allocation param¬ 
eters (e.g., x/, P/ j over X2 interface. If Algorithm 3 executes 
T < Tmax iterations, the overall signalling overhead (e.g., 
number of messages exchanged) is given by 

TV 2 + TV + 2 

nr i = T(0, + l) = Tx- 2 . (52) 

Likewise, for the second case (e.g., when TV > Ui), Algorithm 
2 terminates when there are no unallocated UEs with less 

N (n) 

than their quota requirement (e.g., x„, < n ui ). Hence the 

n — 1 

maximum number of iterations T max = Ui, and the numnber 
of messages exchanged in Algorithm 2 for the second case is 
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given by 


TABLE II 

Summery of performance results 


fli 


■L max 

E ( N - 

2=1 

(N + 1 )Ui 


u, 

i + l) = ^2 N -i + 1 
2=1 

Ui(Ui + 1 ) 

2 


(53) 


Therefore, the overall signalling overhead for the second case 
can be expressed as follows: 


np 08 =Tx 2 ( iV+1 ) [/i U l (U l + 1 ) + 2 (54) 

As can be seen from (l52l > and (l54l) . the overhead of signalling 
increases with the number of UEs and the available RBs. How¬ 
ever, the proposed distributed approach has significantly lower 
computational complexity than the centralized approach (linear 
compared to cubic) and it offers performance improvement 
over the existing solutions (see Section IVIII-Bl i. 


VIII. Performance Evaluation 
A. Simulation Setup 

We develop a discrete-time simulator in MATLAB and 
evaluate the performance of our proposed solution. We sim¬ 
ulate a single three-sectored cell in a rectangular area of 
700 m x 700 m, where the eNB is located in the center of the 
cell and three relays are deployed, i.e., one relay in each sector. 
The CUEs are uniformly distributed within the relay cell. The 
D2D UEs are located according to the clustered distribution 
model |42| . In particular, the D2D transmitters are uniformly 
distributed over a radius D r d\ and the D2D receivers are 
distributed uniformly in the perimeter of the circle with radius 
Dd t d centered at the corresponding D2D transmitter (Fig. 2 in 
|Z|). Both D r d and l),i.d are varied as simulation parameters 
and the values are specified in the corresponding figures. The 
simulation results are averaged over 200 network realizations 
of user locations and channel gains. We consider a snapshot 
model and all the network parameters are assumed to remain 
unchanged during a simulation run. For propagation modeling, 
we consider distance-dependent path-loss, shadow fading, and 
multi-path Rayleigh fading (see Section VII-A in |2]|). 

We measure the uncertainty in channel gains as percent¬ 
ages and assume similar uncertainty bounds in the CQI 
parameters for all the UEs. For example, uncertainty bound 
£ = = £21 = £ 3 ") = £ 4 ”'* = 0.25 refers that uncertainty 

(e.g., estimation error) in the CQI parameters for \/ui,n,l is 
not more than 25% of their nominal values. The simulation 
parameters are similar to those in Table II in {2|. 


B. Results 

A summary of the observations from the performance 
evaluation results is provided in Table HI] 

1) Convergence and goodness of the solution: In Fig. [T] 
we show the convergence behavior of our proposed distributed 
algorithm. In particular, we plot the average achievable data 
rate for the UEs in different network realizations versus the 
number of iterations. The algorithm starts with uniform power 
allocation over RBs, which provides a higher data rate at the 


Observations 

• The proposed distributed solution converges to a 
stable data rate within a few iterations and performs 
close to the optimal data rate with significantly less 
computational complexity. 

• As the distance between D2D peers increases, the 
data rate for direct communication decreases. In 
such cases relaying of D2D traffic can improve the 
end-to-end data rate between D2D peers. 

• After a distance threshold, relay-aided D2D com¬ 
munication provides considerable gain in teims of 
the achievable data rate for the D2D UEs. Even for a 
relatively large distance between the relay node and 
a D2D UE, relaying can provide a better data rate 
compared to direct communication for distant D2D 
peers. There is also a trade-off between achievable 
data rate and robustness against channel uncertainty. 


Ref. figure(s) 
Figs. 1 1121 


Fig. DO 


Figs. 14151 



Fig. 1. Convergence of the proposed solution where the number of CUEs and 
D2D UEs served by each relay node is 5 and 3, receptively (e.g., \Ui | = 8 ). 
D r d and /.),/ 4 are set to 50 m, and uncertainty in CQI parameters is assumed 
to be not more than 25%. 


first iteration; however, it may cause severe interference to 
other receiving nodes. As the algorithm executes, the alloca¬ 
tions of RB and power are updated considering the interference 
threshold and data rate constraints. From this figure it can be 
observed that the solution converges to a stable data rate very 
quickly (e.g., in less than 10 iterations). 

We compare the performance of our proposed scheme with 
a dual-decomposition based suboptimal resource allocation 
scheme proposed in lf43l . We refer to this scheme as existing 
algorithm. In this scheme, the relay node allocates RBs 
considering the data rate requirement and the transmit power 
is updated in an iterative manner by updating the Lagrange 
dual variables. For details refer to |43] Algorithm 2]. The com¬ 
plexity of this algorithm is of O (NUi log N + N log Ui + A), 
where A denotes the number of iterations it takes for the power 
allocation vector to converge l43l . 

In Fig. |2(a)[ we show the performances of the proposed 
distributed scheme and the existing algorithm, and the upper 
bound of the optimal solution which can be obtained in a 
centralized manner using Algorithm [l] We use the MATLAB 
optimization toolbox to obtain this upper bound. We plot 
the average achievable data rate for the UEs versus the 
total number of UEs. The average data rate is given by 
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(a) 


Proposed scheme 
Existing algorithm 


15 18 21 24 27 30 33 

Total number of UEs 


(b) 

Fig. 2. (a) Average achievable data rate for optimal upper bound, distributed 
stable matching and existing algorithm, (b) Efficiency of the proposed solution 
and the existing algorithm. Total number of UEs (i.e., C + D) are varied from 
9 + 6 = 15 to 15 + 18 = 33. D r d and Dd d are assumed to be 50 m. 


E R u h 

R avg = -> w h ere K ch is the achievable data rate 

for UE u. Note that, for a given number of RBs, increasing 
the number of UEs decreases the data rate. Recall that, the 
complexity of both the proposed and reference schemes is 
linear with the number of RBs and UEs; and for the optimal 
solution, the complexity is cubic to the number of RBs and 
UEs. As can be seen from this figure, the proposed approach 
outperforms the existing algorithm and performs close to the 
optimal solution. 

In order to obtain more insights into the performance, in 
Fig |2(b)| we plot the efficiency of the proposed scheme and 
existing algorithm for different number of UEs. Similar to 
1441 Chapter 3], we measure the efficiency as r?(.) = ) , 

where R op tm is the network sum-rate for optimal solution. 
The parameters R prop and R eX st denote the data rate for the 
proposed and existing schemes, respectively, which are used to 
calculate the corresponding efficiency metric r] pmp and rj exs t- 
The closer the value of r]t.) to 1, the nearer the solution is 
to the optimal solution. Clearly, the efficiency of the existing 
algorithm is lower compared to the proposed scheme. From 
the figure we observe that even in a dense network scenario 
(i.e., C + D = 15 + 18 = 33) the proposed scheme performs 


80% close to the optimal solution (compared to 60% for the 
existing algorithm); however, with much less computational 
complexity. 

2 ) Impact of relaying: We compare the performance of the 
proposed method for relay-aided D2D communication with a 
conventional underlay D2D communication scheme. In this 
scheme 0, an RB allocated to CUE can be shared with at 
most one D2D link. The D2D UEs share the same RB(s) (allo¬ 
cated to a CUE using Algorithm [3| and communicate directly 
with their peers without a relay if the data rate requirements 
for both the CUEs and D2D UEs are satisfied; otherwise, the 
D2D UEs refrain from transmitting. We refer to this underlay 
D2D communication scheme 0 as the reference scheme. 
Notice that this conventional (e.g., direct) D2D communication 
approach can save half of the RBs. Therefore, as mentioned in 
footnote |3| data rate of ui in the reference scheme is given by 

R r ff = J2 x^Brb log 2 (1 + Pui ] ful ] \ On the contrary, 

n =1 _ V ' 

data rate of ui in the proposed relay-aided approach is given 
by i?pr°p = £ i x^Brb log 2 (l + l). 

n —1 v 7 

(i) Average achievable data rate vs. distance between D2D 
UEs: The average achievable data rates of D2D UEs for 
both the proposed and reference schemes are illustrated 
in Fig. [3 Although the reference scheme outperforms 
when the distance between the D2D UEs is small (i.e., 
d < 40 m), our proposed approach, which uses relays for 
D2D traffic, can greatly improve the data rate especially 
when the distance increases. This is due to the fact 
that when the distance increases, the performance of 
direct communication deteriorates due to increased signal 
attenuation. Besides, when the D2D UEs share resources 
with only one CUE, the spectrum may not be utilized 
efficiently, and therefore, the achievable rate decreases. 
As a result, the gap between the achievable rate with our 
proposed algorithm and that with the reference scheme 
widens when the distance increases. 

(ii) Gain in aggregate achievable data rate vs. varying dis¬ 
tance between D2D UEs: The gain in terms of aggregate 
achievable data rate under both uncertain and perfect 
CQI is shown in Fig. Q] We calculate the rate gain 
as follows: R gain = Fiplop x 100%, where R plop 
and f? re f denote the aggregate data rate for the D2D 
UEs in the proposed scheme and the reference scheme, 
respectively. The figure shows that, compared to direct 
communication, with the increasing distance between 
D2D UEs, relaying provides considerable gain in terms 
of achievable data rate and hence spectrum utilization. 
As expected, the gain reduces under channel uncertainty 
since the algorithm becomes cautious against channel 
fluctuations and allocates RBs and power accordingly 
to protect the receiving nodes in the network. Note that 
there is a trade-off between performance gain and robust¬ 
ness against channel uncertainty. For example, when the 
distance l) r .d = 50 m, the performance gain of relaying 
under perfect CQI is 30%. In the case of uncertain CQI, 
the gain reduces to 24% and 16% for the uncertainty 
bound parameter £ = 0.25 and £ = 0.50, respectively. 
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Fig. 3. Gain in average achievable data rate with varying distance between 
D2D peers using a setup similar to that for Fig. |T| The reference scheme is 
an underlay D2D communication approach proposed in (5). 



Fig. 4. Gain in aggregate achievable data rate for both perfect and uncertain 
CQI parameters. For uncertain CQI, uncertainty bound £ = 0.25 and 
£ = 0.50 mean that uncertainty in CQI parameters is not more than 25% 
and 50%, respectively. For both the perfect and uncertain cases, there is a 
critical distance beyond which relaying of D2D traffic provides significant 
performance gain. 

As the uncertainty bounds increase, the system becomes 
more roust against uncertainty; however, the achievable 
data rate degrades. 

(iii) Effect of relay-UE distance and distance between D2D 
UEs on rate gain: The performance gain in terms of 
the achievable aggregate data rate under different relay- 
D2D UE distances is shown in Fig. 0 It is clear from 
the figure that, even for relatively large relay-D2D UE 
distances, e.g., D r d > 60 m, relaying D2D traffic 
provides considerable rate gain for distant D2D UEs. 

IX. Conclusion 

We have investigated the radio resource allocation problem 
in a relay-aided D2D network considering uncertainties in 
wireless channels and provided an iterative distributed solution 
to this problem using stable matching. We have analyzed the 
stability, uniqueness, and optimality of the proposed solu¬ 
tion. We have also analyzed the complexity of the proposed 
approach. Numerical results have shown that the distributed 



Fig. 5. Effect of relay distance on rate gain: |C| = 15, 'D = 9. Uncertainty 
in CQI parameters is assumed to be not more than 25%. For every D r ,/, there 
is a distance threshold (i.e., upper position of the light-shaded surface) beyond 
which relaying provides significant gain in terms of aggregate achievable rate. 

solution is close to the centralized optimal solution with sig¬ 
nificantly lower computational complexity. We have compared 
the proposed relay-aided D2D communication scheme with 
an underlay D2D communication scheme. As an extension 
of this work, design and analysis of a unified mdium access 
control (MAC) protocol incorporating mode selection, device 
discovery, and such relay-aided D2D communication in the 
context of LTE-A network will be worth investigating. 

Appendix 

A. Proof of Proposition \2\ 

Note that any arbitrary matching is not necessarily stable. In 
the following, we show that for any given preference profiles, 
each iteration of Algorithm [3] ends up with a stable matching 
(i.e., there is no blocking pair). We prove the proposition by 
contradiction. Let /// be a matching obtained by Algorithm [2] 
at any step t of Algorithm [3] Let us assume that RB n is not 
allocated to UE ui, but it has a higher order in the preference 
list. According to this assumption, the (ui , n ) pair will block 
Pi- 

Since the position of ui in the preference profile of n is 
higher compared to the user in that is matched by pi, i.e., ui A 
pi{n), RB n must select ui before the algorithm terminates. 
However, the pair ( ui, n) does not match each other in the 
matching outcome pi. This implies that iq rejects n (e.g., line 
7 in Algorithm [2j and ( iii,n ) is a better assignment. As a 
result, the pair ( ui , n) will not block pi, which contradicts our 
assumption. Consequently, the matching outcome pi leads to 
a stable matching since no blocking pair exists and the proof 
concludes. 

B. Proof of Proposition [7] 

The proof is followed by the induction of number of 
users Ui, that are supported by relay l. For instance, let 
k ui = 1, Vui € Ui (the proof for k ui > 1 can be done 
analogously introducing dummy rows [i.e., UEs] in the utility 
matrix). The basis (i.e., Ui = 1) is trivial, since the only user 
definitely gets the best RBs according to her preference. When 
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Ui > 2, let us consider R\ to be the maximal entity of the 
utility matrix 11;. For instance, let the matrix 11; be obtained 
by removing the z-th row and j-th column from the utility 
matrix if;. If Hi is a stable matching for if;, then by definition 
Hi (i) = j and hence hi \ {(/, j)} must be a stable matching 
for if;. By induction, there exists a unique stable matching hi 
for the smaller matrix if;. Therefore, the proof is concluded 
due to the fact that /i; = //; U {(*,_))} is the unique stable 
matching for the utility matrix if;. 


C. Proof of Proposition [7] 

Without loss of generality, let 93 Ui! ;(/i;) denote the data rate 
achieved by UE it; for any matching /z; for given uncertainty 
bounds and 93; (/i;) = 93 Uij ;(/i;) is the sum-rate of all 

ui£Ui 

UEs. On the contrary, let /i/ denote an arbitrary unstable 
outcome better than /q, i.e., pi can achieve a better sum- 
rate. There are two cases that make /}; unstable: 1) lack of 
individual rationality, and/or 2) blocked by a UE-RB pair 0. 
We analyze both the cases below. 

Case 1 (lack of individual rationality): If RB n is not 
individually rational, then the utility of n can be improved 
by removing user /2;(n) with any arbitrarily user it; = pi(n). 
Hence, the utility of m increases and 93 Uli ;(/i;) < 93 Uij ;(/i;). 

Case 2 (pi is blocked): When /T; is blocked by any UE-RB 
pair (it;,n), RB n strictly prefers UE it; to pi(n) and one of 
the following conditions must be true: 

(i) ui strictly prefers n to some n G jtt;(it;), or 

(ii) p,i(ui) < tt Ul and n is acceptable to it;. 

If condition (i) is true, we can obtain a stable matching /z; 
by interchanging n and n for it; as follows: 

m(ui) = {jh(ui)\n}Un. (A.l) 

Hence, the new data rate of UE it; is 

KuM= E R W = R % ) + E 

>R% ) + E R u}= E R $=n*M 

(A.2) 

where Ru' 1 is given by (l45l >. Since it; strictly prefers RB n 
to n and the data rates for other UEs remain unchanged, for 
condition (i), it can be shown that 93;(/i;) > 93; (//;). 

When condition (ii) is true, 

9W(w)= E R W+ R W 

jePi(ui) 

> E =K uu i(jh). (A.3) 


Let iii = ft; (it) with data rate R7\ Then 

KuM= E 

iefipui) 

< E (A.4) 

jeih(ui) 

From ( IA.3b and (1A.4I >. neither 93; (/i;) > 93;(/£;) nor 93; (hi) > 
93 i(pi). Since for both cases 1 ) and 2) there is no outcome 
Hi better than //; , by Definition [6] hi is an optimal allocation 
and the proof follows. 
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